AD-A163  825 
UNCLASSIFIED 


MODIFICATION  OF  THE  IONOSPHERIC  ELECTRON  VELOCITY 
DISTRIBUTION  FUNCTION  D  (U)  CALIFORNIA  UNIV  LOS 
ANGELES  CENTER  FOR  PLASMA  PHYSICS  AND  FUS 
M  SHOUCRI  ET  AL  DEC  85  PPG-913 


1/1 

ML 


F/G  4/1 


MICROCOPY  RESOLUTION  TEST  CHART 

NATjOWAL  BUREAU  OF  STAN0AR0S-I963-A 


7TS  T7S  ~  W  TS^  -r^ir^TTWV’lTEVKlCVTI  1 


MODIFICATION  OF  THE  IONOSPHERIC 


ELECTRON  VELOCITY  DISTRIBUTION  FUNCTION 


DUE  TO  RESONANT  ABSORPTION  OF  HF  WAVES 


Merit  Shoucri,  G.J.  Morales,  &  J.E.  Maggs 


December  1985 


PPG-913 


CENTER  FOR 

PLASMA  PHYSICS 
AND 

FUSION  ENGINEERING 


UNIVERSITY  OF  CALIFORNIA 
LOS  ANGELES 


DTIC 

ELECTE 

JAN  8  1986 


DBTMBUTIOW  STATEMENT  A 


Approved  tot  public  release) 
DUtdbutioo  Unlimited  v'“ 


MODIFICATION  OF  THE  IONOSPHERIC 


ELECTRON  VELOCITY  DISTRIBUTION  FUNCTION 
DUE  TO  RESONANT  ABSORPTION  OF  HF  WAVES 


Merit  Shoucri,  G.J.  Morales,  &  J.E.  Maggs 


December  1985  PPG-913 

Contract  N00014-7 5-C-0476 


University  of  California 
Department  of  Physics 
Los  Angeles,  CA  90024 


I 

s 


DTIC 

ELECTE 
JAN  8  1986 


B 


nPaffMBUTlOW  STATEMENT  A~\ 

Approved  toe  public  release!  f 
Distribution  Unlimited  jyj 


MODIFICATION  OF  THE  IONOSPHERIC  ELECTRON  VELOCITY 


DISTRIBUTION  FUNCTION  DUE  TO  RESONANT  ABSORPTION  OF  HF  WAVES 


Merit  Shoucri,  G.J.  Morales,  and  J.E.  Maggs 
Physics  Department 

University  of  California  at  Los  Angeles 
Los  Angeles,  CA  90024 


ABSTRACT 

\ 

The  nonlinear  modification  of  the  electron  distribution  function  caused 
by  the  interaction  with  a  resonant  electrostatic  field  is  studied  analytically 
using  a  perturbation  analysis.  This  interaction  produces  tail  heating,  and 
for  a  range  of  parameter  values,  can  form  a  bump  in  the  distribution,  which 
leads  to  emission  of  secondary  waves.  These  results  are  applied  to  the 
ionosphere's  F  and  upper  E  regions  in  the  context  of  ionospheric  modification 
experiments  using  ground-based  powerful  HF  transmitters.  Analogies  becween 
our  results  and  experimental  measurements  are  discussed.  'J-y  J  I 
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INTRODUCTION 


The  interaction  of  powerful  high  frequency  (HF)  electromagnetic  waves 
with  ionospheric  plasmas  can  lead  to  the  generation  of  electrostatic  plasma 
waves  through  various  processes,  e.g.,  linear  mode  conversion,  direct 
conversion,  parametric  instabilities,  etc.  These  processes  occur  at  a  height 
where  the  angular  frequency  of  the  HF  wave  u0  nearly  matches  the  local 
electron  plasma  frequency  oip.  The  relatively  slow  phase  velocity  of  the 
excited  electrostatic  modes  allows  the  possibility  of  strong  interactions  with 
ionospheric  electrons.  These  interactions  can  result  in  significant  modifi¬ 
cations  of  the  background  electron  distribution  function.  The  signatures 
of  such  Interactions  are  often  observed  during  ionospheric  modification 
experiments  [see,  for  reference,  volume  j>,  issue  11,  of  Radio  Science  (1974), 
and  volume  44,  issue  12,  of  Journal  of  Atmospheric  and  Terrestrial  Physics 

(1982)].  For  example,  during  HF  heating  experiments,  fast  electrons  produce 
enhanced  airglow  at  6300A  and  5577A,  caused  by  the  excitation  of  atomic  oxygen 

lines  [Sipler  and  Blondi,  1972;  Haslett  and  Megill,  1974;  Carleson,  et  al., 

1982].  Non-ohmic  heating  of  the  ionosphere  which  is  partially  attributed  to 

tail  particles  produced  during  the  UF^heating  has  also  been  observed  [Meltz, 

et  al.,  1974;  Mantas,  et  al.,  1981].  The  purpose  of  the  present  work  is  to 

calculate  the  modification  in  the  electron  velocity  distribution  function 

resulting  from  the  resonant  excitation  of  electrostatic  oscillations  in  the 

ionosphere.  It  is  hoped  that  the  analytic  results  may  help  illuminate  some 

of  the  various  experimentally  observed  phenomena,  and  perhaps  suggest  new 

experiments  which  would  probe  the  nature  of  the  electron  acceleration 


mechanisms 
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This  study  does  not  address  the  different  mechanisms  that  can  generate 
plasma  waves.  Such  studies  are  available  In  the  literature.  For  example,  the 
generalized  mode  conversion  problem  Is  treated  In  detail  by  Mj^lhus  [1983]; 
direct  conversion  Is  analyzed  by  Wong,  et  al.,  [1981];  and  parametric 
Instabilities  In  the  Ionosphere  are  reviewed  by  Fejer  [1979].  The  present 
work  assumes  a  given  spatial  structure  for  the  electric  field  near  the  plasma 
resonance,  and  proceeds  to  calculate  the  resultant  kinetic  Interactions.  Two 
typical  resonant  structures  are  considered:  one  corresponds  to  the  convective 
limitation  of  the  plasma  resonance  (Airy  pattern),  and  the  other  Is 
appropriate  to  a  colllslonally  limited  resonance  (Lorentzian  pattern). 

Figure  1  Illustrates  the  geometry  of  the  model  problem,  and  shows  the 
relevant  scales.  The  localized  electrostatic  field  Ej(z)  is  considered  to  be 
parallel  to  the  earth's  magnetic  field  if0  at  z  *  0.  Its  characteristic  width 
Is  Ax*  ~  LnCkoLn)-273  for  the  Airy  pattern,  or  Az*  ~  (rLn)  for  the  Lorentzian 
pattern,  where  T  -  ve/o>0,  ve  is  the  electron  collision  frequency,  Ln  the 
density  scalelength  along  B0,  and  kp  is  the  Debye  wavenumber.  Electrons 
Impinging  upon  the  localized  field  from  the  high  density  side  can  be 
accelerated  as  they  pass  through  the  interaction  region  Az^,  hence,  modifying 
the  distribution  function  downstream  at  z  ■  z0.  This  modified  distribution 
changes  slowly  within  a  slowing  down  length  L(u)  for  those  electrons  with 
velocity  u  parallel  to  the  magnetic  field;  beyond  L(u),  the  distribution 
relaxes  to  an  appropriate  equilibrium. 

In  this  work,  we  analytically  calculate  the  modified  electron  distribu¬ 
tion  existing  in  the  range  between  z0  and  L(u)  by  using  a  perturbation 
analysis  correct  to  second  order  In  the  strength  of  the  electric  field.  The 


general  result  we  derive  Is  applied  to  a  specific  bi-Maxwellian  background 
electron  distribution,  characteristic  of  an  ionospheric  F  region.  Parameters 
of  the  background  plasma  are  varied  to  assess  the  effect  of  different 
ionospheric  conditions. 

An  interesting  feature  found  in  the  analysis  of  the  Airy  and  Lorentzian 
patterns  is  that  both  structures  have  a  broadband,  unidirectional  Fourier 
spectrum,  with  identical  power  spectra  characterized  by  a  colllslonally 
dependent  exponential  envelope.  For  electric  field  amplitudes  above  a  certain 
threshold,  a  bump  can  be  formed  in  the  tall  of  the  velocity  distribution 
function  downstream  from  the  resonance  layer.  The  region  of  positive  slope  in 
the  velocity  distribution  may  trigger  the  growth  of  a  broad  spectrum  of 
sideband  waves  in  the  neighborhood  of  <%.  The  width  of  the  sideband  frequency 
spectrum  is  determined  by  the  electric  field  strength,  its  effective  phase 
velocity,  and  by  the  background  plasma  parameters.  It  is  worth  noting  that 
these  sidebands  are  analogous  to  ones  observed  in  the  laboratory  by  Starke  and 
Malmberg  [1978].  Another  interesting  result  obtained  from  the  calculation  is 
the  observation  of  heating  of  particles  with  velocities  of  the  order  of,  or 
larger  than,  the  characteristic  phase  velocity  of  the  resonant  field.  These 
tall  electrons  ultimately  loose  their  energy  to  neutral  particles  and  bulk 
electrons  by  different  collisional  processes  that  occur  in  the  ionosphere. 

The  observed  airglow  enhancement  is  but  one  of  these  processes.  Through  them 
one  can  estimate  the  slowing  down  length  for  fast  electrons. 

The  description  and  analysis  of  the  model  are  presented  in  section  2. 
Results  obtained  from  the  model  are  discussed  in  section  3,  while  conclusions 
are  given  in  section  4. 
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2.  MODEL 


This  analysis  considers  a  nonuniform,  weakly  collisional  (ye  «  «-b), 
magnetized  plasma  characteristic  of  the  ionospheric  F  and  upper  E  regions. 
The  resonantly  excited  electric  fields  are  assumed  to  be  primarily  aligned 

A 

along  the  magnetic  field  £B0,  i.e.. 


E(r,t)  a  z,  [E|(z)exp(-iu0t)  +  c.c.]  ; 


(1) 


hence,  the  Boltzmann  equation  is  solved  by  separating  the  total  distribution 
function  f(r,v,t)  into  an  unperturbed  symmetric  perpendicular  distribution 
F0(vj.)  and  a  perturbed  parallel  (field  aligned)  distribution  g(z,u,t). 


f(r,v,t)  -  FoUjJgCz.u.t), 


(2) 


^  ^  ^ 
where  v  *  vj^  +  zu,  and  |  v^|  “  v^. 


2.1  Perturbation  Analysis 


We  proceed  with  a  perturbation  calculation  correct  to  order  | E | 2 . 
Equation  (2)  is  therefore  expanded  according  to 


f(r,v,t)  -  fQ(r,v)  +  fj(r,v,t)  +  f2(r,v,t)  +  .  .  .  . 


(3) 


The  spatially  dependent  equilibrium  distribution  is 


Mr.*)  "  F  (v ,  )g  (u)n  (z) , 


nQ(z)  -  S0(l  -  z/Lfl), 


(4) 


where  a  linear  profile  is  assumed  around  the  region  of  interest,  and  N0  is  the 
density  at  the  plasma  resonance  (z  ■  0).  The  normalization  conditions  for 


f- f.  ■r.-r.  V- 


the  zero  order  distribution  are  /du  g0(u)  ■  1  and  /d2vjF0(vp  -  1.  Symmetry 
in  the  perpendicular  direction  also  implies  that  /d2vjvjF0( v^)  “  0»  The  first 
order  perturbed  distribution  oscillating  at  frequency  (Uq  is 


f^r.v.t)  -  Fo(vl> lg1(*»u)exp(-ia»ot)  +  c.c.] 

and  represents  Che  linear  response  of  the  plasma.  The  second  order  perturbed 
distribution  function  contains  both  a  time  independent  term,  corresponding  to 
the  modification  of  the  equilibrium  by  the  resonant  field,  and  a  second 
harmonic  response 

f2(r,v,t)  -  F^Vj)  {5g2(z,u)  +  [g2(z,u)exp(-i2unjt)  +  c.c.]}. 

However,  for  the  purpose  of  the  present  study,  the  second  harmonic  term 
will  not  be  considered  and  is  subsequently  neglected. 


The  Boltzmann  equation  is  written  as 


_3_ 

3z 


n(z)exp(-iuot)  + 


,xzB 

l—o 


].V+}  Fo(Vi)g(z,u,t)  - 


(if) 


3t  Jcoll* 


(5) 


The  unperturbed  equilibrium  is  determined ' by 


3n  (z) 

u  _2 - g0(u)Fo(vi) 


(— ) 

'•St J 


coll,o 


(6) 


which  states  that  the  plasma  density  gradient  is  determined  by  a  balance 
between  the  sources  and  the  transport  mechanisms,  both  plasma  and  atomic, 
existing  in  the  ionosphere  (l.e.,  elastic  and  inelastic  collisions,  ionlza- 
tlon,  and  recombination  processes)  and  contained  in  the  term  (-s^)co^  Q. 
Since  the  present  study  aims  to  resolve  the  acceleration  mechanisms  in  the 
neighborhood  of  plasma  resonance,  a  linear  approximation  (4)  for  the 


unperturbed  density  profile  is  made,  with  the  scalelength  given  by 


L 

n 


3no ,-1 

3z  J 


(7) 


The  first  order  parallel  distribution  function  is  obtained  from  (5)  by 
integrating  over  perpendicular  velocities 


3gj(z,u) 

3z 


-  i<g1(z,u) 


en  (z)E  (z)  3g  (u) 

O  8  Q 

mu  3u  * 


(8) 


u>0  +  iv(u) 

where  <  - - ,  and  the  first  order  linearized  contribution  of  the 

u 

right-hand  side  term  in  (5)  is  approximated  by  -v(u)g  (z,u), 
v(u)  -  fd  vj^v(v)F0(vj_)  is  the  slowing  down  rate  for  electrons  with  velocity 
u  averaged  over  perpendicular  velocities  (not  to  be  confused  with  the  average 
collision  frequency  for  bulk  electrons  ve).  Equation  (8)  is  a  linear 
differential  equation  in  z  which  is  integrated  along  its  characteristics. 
Separating  right-going  [g|  (z,u)]  from  left-going  [g1  (z,u)]  particles: 


S 


i 

1 


(z,u,t) 


e  *o  (u)  r 
mu  3u 


dk  E(k)exp[i(kz  -  a>ot)] 
2ffi  (k  -  ic) 


Mo  i 
(k  -  K)LnJ  ’ 


(9) 


where  E(k)  is  the  Fourier  transform  of  the  field  Ej(z),  and  in  general 
gn(u)  ■  g^(u)9(u)  +  gn(u)9(-u),  8(u)  being  the  Heaviside  unit  step  function. 


Carrying  out  the  analysis  to  second  order  in  the  resonant  field,  the  time 
averaged  component  of  (5)  integrated  over  perpendicular  velocities  is  written 


as: 


u  Iz  5g2(z,u) 


3g,(z,u)  *  3g  (z,u) 

a  [£,<z)  - K - +  - 3^ - ] 


3u 


9u 


(10) 


/d2vi<f'f)coll,2>’ 


where  the  time-averaged  second  order  distribution  is  <f2(z,v, t)>  * 

3f 

6g2(z,u)Fo(v^),  and  2^  rePresenta  the  time-averaged  second  order 

collisional  processes.  The  interaction  between  the  particles  and  the  resonant 
electric  field  Eg(z)  occurs  within  a  spatially  localized  region  Az^  (see 
Figure  1)  that  is  assumed  to  be  much  smaller  than  the  slowing  down  length 
L(u)  *  u/v(u)  for  velocities  of  interest.  In  the  ionosphere's  F  region  and 
the  upper  E  region,  this  assumption  is  well  justified.  However,  for  lower 
altitudes,  the  stopping  distance  decreases  with  increasing  neutral  density  and 
this  assumption  can  breakdown.  It  is  also  understood  that  the  density 
scalelevtgth  is  also  large,  hence  we  write 


Ln,  L(u)  »  jAzJ  .  (11) 

For  distances  from  the  resonant  layer  z  <  L(u),  Sg2(z,u)  can  be  obtained  by 
neglecting  the  effect  of  collisions,  i.e.,  the  right-hand  side  of  (10),  and  we 
can  write: 


<5g~  ( z ,  u) 


^_(z)  3 


16ff3mu  u 


.  3g~(u)  » 

—  f-  — -  / 

a.,  1..  3u  J 


dk 


E»(k) 

(k-K*) 


/  dk’ 


E(k’) 

(k’-k) 


exp 


[iz(k  -  k*) 


iN  /n  (z)  IN  /n  (z) 
r  1  j.  0  0  O  O  1  .  1 

(k-**)L  ■  (k’-k)L  J  c*c. 

n  n 


(12) 


In  deriving  (12),  we  have  made  the  approximation  that  L(u)  +  «  for  5g2±. 
This  approximation  is  in  accordance  with  (11),  which  expresses  the  localized 
nature  of  Ej(z).  Furthermore,  since  Ln  >>  |AzjJ,  the  last  two  terms  in  (12) 
make  a  small  contribution,  and  can  thus  be  neglected. 


Up  to  this  point,  no  specific  form  has  been  assumed  for  either  g~(u) 

or  E(k)  except,  of  course,  that  they  must  be  physical  (i.e.,  that  they  are 
well  behaved  at  infinity).  More  specifically,  we  have  stated  that  Ej(z)  is  a 
spatially  localized  field  near  the  plasma  resonance  region,  which  implies  that 
its  Fourier  spectrum  E(k)  is  broadband.  This  being  the  case,  if  we  further 
require  that  the  spectral  amplitude  | E(k) |  of  the  electric  field  is  an  entire 
function  (i.e.,  does  not  have  singularities  in  the  finite  k-plane)  then,  for 
this  second  order  analysis,  one  can  average  the  effect  of  the  phases  by 
calculating  5g^(z,u)  in  the  asymptotic  region  where  z  >_  zQ  »  Az*  (see  Figure 
1).  Equation  (12)  reduces  to  a  single  integral 


<5g~  (z  ,u) 


.  »  [»  3g»  (U)  f*  IE-OQ] 2  ,  . 

16ir^mu  3u  *-u  3u  *  i(k-ic)  *  ‘ 


(13) 


Using  (12)  or  (13),  the  total  time-averaged  electron  distribution  correct 
to  second  order  is: 

<  f  (z , v, t)  >  -  [nQ(z)g~(u)  +  5g'(z,u)j  Fq  (Vj).  (14) 

2.2  Electrostatic  Field  Structure 


To  give  a  concrete  example  of  the  changes  that  occur  in  the  electron 
distribution,  we  consider  two  field  patterns,  the  Airy  function  and  the 
Lorentzlan.  These  two  field  structures  are  characteristic  of  resonantly 


excited  plasmas.  The  Airy  structure  arises  when  linear  mode  conversion  limits 
the  resonance  amplitude  excited  by  an  incident  electromagnetic  wave  propaga¬ 
ting  into  a  plasma  with  a  density  gradient.  The  Lorentzian  pattern  replaces 
the  Airy  pattern  whenever  the  mean  free  path  of  the  thermal  electrons  is  smal¬ 
ler  than  the  characteristic  Airy  scalelength,  i.e.,  (v  /v  )  <  L  (k_L  )“2/3,  v 
being  the  electron  thermal  velocity. 


The  electric  field  for  an  electrostatic  plasma  wave,  propagating  along 
the  magnetic  field  in  the  positive  z  direction  is  given  by  [Morales  and  Lee, 
1974] 


2/3 


E„(z)  -  -iirEQ(kDLn)  [Ai(-5  -  ir')  -  iGi (-C  -  iP)]  , 


(15) 


where  C  *  (k^L^)2''3  z/Ln,  P  ■  (kQLn)2/3r,  T  »  v^/ tn^  Ai  and  Gi  are  solutions 
to  Airy's  equation  given  in  Abramowitz  and  Stegun  [1965],  and  E0  is  a  field 
amplitude  that  is  related  by  a  mode  conversion  efficiency  to  the  incident 
electromagnetic  wave  power.  In  the  limit  where  collisions  dominate,  we  obtain 
the  Lorentzian  profile 


E,U) 


E  L 
o  n 


z  +  ian  * 


(16) 


It  is  worth  noting  that  for  z/Ln  r,  (15)  is  approximated  to  leading  order  by 
(16).  This  effect  is  illustrated  In  Figure  2a,  which  depicts  the  magnitudes 
of  both  the  Airy  and  Lorentzian  patterns  for  two  values  of  the  damping  factor, 
P  *  0.2  and  0.8.  Figure  2b  displays  the  real  parts  of  both  patterns  for 
P  *  0.4.  The  Lorentzian  consists  mainly  of  a  single  wavelength,  that  is 
almost  in  phase  with  the  Airy  pattern  near  the  origin,  and  travelling  down  the 
density  gradient.  This  feature  is  also  found  by  taking  the  Fourier  transform 
of  both  field  structures: 


-  Sl.w.it..**...* I-*..*,  ,"L  X  -m. 


Wy 


a.  A  a  M  - 


+  }]} . 


(17) 


E  (k) 

A 


2irE  L  0(k)exp 
o  n 


<kLn> 

3(Vn 


Vk)  -  2irEoLn8(k)exp  [-  -  |i]. 


(18) 


where  again  9(k)  is  Che  Heaviside  unit  step  function,  E^(k)  is  Che  Fourier 
transform  of  the  Airy  pattern  of  (15)  and  E^(k)  that  of  the  Lorentzian  of 
(16);  the  characteristic  wavenumber  k0  ■  %/a  *  (rLjj)”1,  where  a  is  the 
effective  phase  velocity  of  the  Lorentzian  and/or  the  Airy  envelope.  The 
derivations  of  (17)  and  (18)  are  shown  in  the  Appendix. 

The  Fourier  transforms  of  (17)  and  (18)  have  common  features.  They  are 
both  broadband  and  unidirectional.  In  fact,  they  only  differ  by  a  highly 
oscillatory  phase  term  in  the  Airy  spectrum.  The  power  spectra  for  the  Airy 
and  the  Lorentzian  are  therefore  identical.  Accordingly,  in  the  limit  of 
second  order  analysis,  a  particle  passing  through  either  field  pattern  emerges 
downstream  at  zQ»  Az^,  away  from  the  interaction  region  having,  on  the 
average,  experienced  identical  effects.  The  magnitude  of  the  electric  field 
Fourier  spectrum  is  plotted  in  Figure  3  as  a  function  of  both  wavenumber  k  and 
parallel  velocity  u,  where  k  ■  u>0/u.  The  resonant  field  k-components  interact 
with  the  background  electrons  and  accelerate  them.  Components  varying  from 
zero  up  to  values  on  the  order  of  kg  -  (H,n)_1  <  kj)  efficiently  interact  with 
the  background  particles.  As  k  increases,  the  interaction  with  the  resonant 
electric  field  structure  decreases  until  k  -*■  kjj,  where  Landau  damping  destroys 
the  wave  completely.  In  terras  of  the  velocity  space  interaction,  this 
translates  as  acceleration  of  particles  mainly  in  the  tail  by  an  electric 
field  of  constant  magnitude  produced  by  the  plateau  that  exists  in  the 
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exponential  behavior  of  the  spectrum  from  ve  <  u  «  a  up  to  u  +  +•.  For 
velocities  u  <  a,  closer  to  the  bulk  of  the  velocity  distribution,  the  field 
has  negligible  k-components  and  particle  acceleration  is  small.  As  u  -*■  v 
Landau  damping  completely  eliminates  any  k-components. 


Substituting  for  either  of  the  field  patterns  in  (13)  gives 


.  +,  .  “p(z)(EoLn>  3  fl  3go  -2a/u  r  /2v(u) . 

5g2(z’u)  *  2SH —  m  t  lu-  *  ^C0S 


+  Ei  (2a/u)sin  )] } 


(19) 


for  u  >  0,  and 
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for  u  <  0;  Ei  and  E*  are  exponential  integral  functions  defined  in  Abramowitz 
and  Stegun  [1965].  For  the  case  v(u)  -*■  0,  i.e.,  negligible  slowing  down  of 
fast  particles,  then  6g^  ♦  0  and 


,  xu  (z)(E  L  )  a  ,  3g  .  , 

■+/  v  m  p  _  o  n  d  |4  °o  e-2a/u 


6g2(z,u) 


2mu 


3u  Lu  3u 


]  • 


(21) 


No  wave  particle  resonances  occur  for  particles  moving  up  the  geomagnetic 
field  (i.e.,  in  the  negative  z  direction)  due  to  the  unidirectional  nature  of 
the  resonant  field  wavenumber  spectrum.  If  the  resonant  field  spectrum  had  k 
components  in  the  direction  opposite  to  the  geomagnetic  field  (negative  z 
direction),  then  particles  with  u  <  0  would  be  accelerated  yielding  a  non  zero 

5g2* 


2.3  Zero  Order  Distribution 


A  self-consistent  calculation  of  the  zero  order  distribution  function 
fQ(z,v)  requires  the  solution  of  equation  (6).  We  shall  not  attempt  this, 
however,  as  it  is  quite  difficult  for  F  region  conditions  and  contributes 
little  to  the  features  we  wish  to  emphasize.  Instead,  we  consider  a  model 
characteristic  of  an  equilibrium  ionosphere  that  can  simulate  different 
regimes  of  Che  ionosphere  by  changing  the  value  of  the  parameters.  Through 
this  approach,  we  exhibit  the  main  features  contained  in  the  present 
analysis. 

A  bl-Maxwellian  models  the  zero  order  electron  distribution  and  is 
expressed  as 
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where  A  is  the  fractional  density  of  the  energetic  tall,  usually  of  the  order 
of  a  few  percent,  ve  is  the  thermal  velocity  of  the  bulk  plasma,  and  vf  is  the 
characteristic  velocity  of  the  tail.  A,  ve,  and  Vf  are  free  parameters  that 
model  the  various  ionospheric  conditions.  Using  (21)  and  (22),  the  resulting 
total  electron  distribution  function  of  (14)  (accurate  to  second  order)  is 
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where  Eg  -  (j  mvf  )/(/2ireLn)  is  a  normalizing  field  that  scales  as  the  average 
tail  energy  per  unit  charge  over  Che  density  scalelength.  If  there  is  no 

Initial  tail,  l.e.,  only  a  single  Maxwellian  background,  then  A  *  1  and 


vf  ■  ve.  Figure  4  shows  a  plot  of  (23)  for  &  *  0.01,  vf  «  10ve,  and  a 
resonant  field  having  an  effective  phase  velocity  a  -  20ve.  For  E0  -  3ES> 
particles  diffuse  toward  higher  velocities  causing  a  swelling  of  the  tail. 

This  tall  heating  is  produced  by  the  power  spectrum  of  the  applied  electric 
field  E(k  ■  (Dq/u)  displayed  in  Figure  3.  The  heating  is  most  significant  for 
large  velocities  where  the  spectral  density  is  largest  and  nearly  constant. 

Equation  (23)  shows  that  the  exponentially  decaying  Gaussian  tail  of  the 

background  is  enhanced  and  spread  out  by  the  perturbation  term.  The  larger 

the  amplitude  of  the  applied  field,  the  larger  the  diffusion  toward  the  tail. 

However,  since  there  are  more  particles  at  lower  velocities,  one  can  see  that 

because  the  magnitude  of  E(k)  Increases  with  velocity  (Figure  3),  there  may  be 

a  velocity  bin  where  more  particles  diffuse  into  it  from  the  bulk  (and  more 

populated)  side  than  there  are  particles  leaving  that  region  toward  higher 

velocities.  This  effect  may  produce  a  localized  swelling  of  the  tail  that 

would  ultimately  cause  a  bump  in  the  distribution  for  large  enough  fields. 

The  region  where  this  process  would  occur  is,  of  course,  a  function  of  both  Vf 

and  a,  i.e.,  the  tail  and  Che  resonant  field's  effective  phase  velocities, 

respectively.  The  condition  for  bump  production  is  analyzed  in  the  following 

section. 

i 

\ 

\ 

2.4  Threshold  for  Bump  Production 

We  now  investigate  the  condition  for  producing  a  region  of  positive  slope 
in  the  distribution  function  of  (23).  For  >U~  2.  one  ®ust  satisfy  the 


condition 
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Again,  for  the  case  of  a  single  Maxwellian  plasma,  a  similar  expression  is 


obtained  with  Vf  *  ve.  Positive  slopes  will  occur  on  the  distribution  only  in 


the  region  where  the  denominator  of  (24)  is  positive,  i.e.,  whenever 
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It  should  be  emphasized  that  although  the  positive  slope  condition  of 


(24)  is  calculated  for  the  specific  spectra  of  the  Airy  and  the  Lorentzian 


field  structures,  this  method  is  applicable  in  general  to  any  spectrum 


satisfying  the  conditions  leading  to  (13).  In  this  case,  the  condition  of 


(24)  for  a  positive  slope  in  the  tail  distribution  is  generalized  to 
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A  is  the  amplitude  of  the  applied  field,  and  w  its  characteristic  width. 


can  be  shown  that  for  a  Lorentzian  or  an  Airy  pattern  (26)  reduces  to  (24) 


2 . 5  Sideband  Emissions 


2.5.1  Threshold  for  Growth 

The  existence  of  a  positive  slope  in  the  distribution  suggests  that 
plasma  waves  with  phase  velocities  resonant  with  particles  in  the  positive 
slope  region  can  become  unstable.  These  unstable  plasma  waves  form  broad 
sidebands  in  the  neighborhood  of  the  main  frequency  ojq.  For  sideband 
amplification,  the  field  E0  must  be  large  enough  not  only  to  produce  a 
positive  slope  in  the  distribution,  but  also  to  overcome  collisional  and 
Landau  damping.  This  condition  is  written  as 
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(27) 


where  the  sideband  frequency  is  ois  *  Mj.  +  iy,  and  we  assume  Y  «  ojj..  Equation 
(27)  is  accurate  to  order  (vg/u  ).  From  (23),  the  condition  of  (27)  is 
rewritten  as: 
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where  again  in  the  case  of  a  single  Maxwellian  plasma  A  ♦  1  and  v^  «  Vg. 
Comparing  (28)  to  (24),  one  concludes  that  the  addition  of  the  term 
proportional  to  T  in  the  former  increases  the  wave  growth  threshold  over  the 
threshold  for  positive  slope.  In  practice,  however,  the  terra  containing  T  in 
(28)  is  usually  small  compared  to  unity  in  the  square  brackets. 
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Given  a  certain  field  E0  above  the  instability  threshold  of  (28),  we 
expect  growth  of  sideband  waves  with  growth  rates  given  by  (27).  The  growth 
rate  is  largest  for  the  sideband  whose  phase  velocity  matches  the  largest 
positive  slope,  and  decreases  away  from  this  value  for  both  larger  and  smaller 
velocities  up  to  frequencies  with  phase  velocities  which  correspond  to 
negative  slope  on  the  distribution. 

2.5.2.  Sideband  Frequency  Spectrum 

The  frequencies  of  the  growing  sidebands  can  be  determined  in  the  case  of 
a  homogeneous  medium,  from  the  dispersion  relation.  For  phase  velocities 
larger  than  the  background  thermal  velocity,  the  dispersion  relation  for 
Bohor-Gross  waves  is  approximately 
2  2  _2  2 

»s  ~  -  3ve/u  ),  (29) 

where  it  is  assumed  that  u»s  is  real.  Velocities  in  the  positive  slope  region 
usually  satisfy  the  condition  u  £  a,  while  the  Airy  electric  field  structure 
at  plasma  resonance  satisfies  the  relation 
2  2  _2  2 

oi  ~  o)  /(I  -  3v  /a  ).  (30) 

o  —  p  e 

With  u  1  a,  equations  (29)  and  (30)  indicate  that  sidebands  should  be  radiated 
at  frequencies  wg  £  o»0,  i.e.,  larger  than  the  initial  resonant  frequency. 

Note  that  for  the  Lorentzian  (cold  plasma),  the  phase  velocity  is  infinite  at 
resonance,  so  that  the  preceeding  statement  is  still  correct.  The  wavenumbers 
corresponding  to  these  waves  are  obtained  by  the  condition  ks  *  u>s/u. 
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For  Che  case  of  an  inhomogeneous  plasma,  recall  Chat  the  distribution 
given  in  (23)  is  valid  in  the  spatial  range  z0  £  z  <  L(u),  where  z0  is  an 
arbitrarily  chosen  location  satisfying  the  condition  z0  >>  Az^  (Figure  1). 
Equation  (29)  is  therefore  rewritten  using  the  position-dependent  plasma 
frequency: 

2  2 

u»  (z)  -  (a  (o)  (  1  -  z/Ln),  (31) 

P  P 

where  <i»p(o)  ■  utp.  Therefore,  in  the  inhomogeneous  plasma,  although  the 
sideband  waves  have  the  same  phase  velocities  as  in  the  homogeneous  case, 
i.e.,  in  the  range  u  1  a,  their  frequencies  are  smaller  because  a>p(z)  <  ojp(o). 
As  a  result,  the  spectrum  of  the  sidebands  is  shifted  to  lower  frequencies 
than  for  a  comparable  homogeneous  case,  as  is  clear  from  (29). 

The  sideband  frequency  spectrum  is  therefore  determined  by  two  factors, 
namely,  that  the  bump  occurs  at  velocities  u  1  a,  and  the  plasma  frequency 
decreases  away  from  the  resonance  region.  Furthermore,  since  the  perturbed 
distribution  is  valid  at  any  z  satisfying  the  condition  zQ  z  <  L(u),  side¬ 
bands  should  be  amplified  for  all  z's  in  this  range  at  frequencies  satisfying 
the  local  dispersion  relation.  Backscattered  Thomson  radar  returns  would 
therefore  see  sidebands  caused  by  the  integrated  effect  from  the  different 
heights.  The  recorded  frequency  spectrum  would  manifest  itself  as  a  bump-like 
distribution  of  broadband  frequencies  that  spread  around  the  initial  frequency 
<u0,  depending  on  the  prevailing  ionospheric  and  resonant  field  parameters 
(ve,  Vf,  a,  zQ,  l^j,  .  .  •). 

We  have  given  a  simple  argument  to  the  effect  of  the  inhomogeneity  on  the 
generated  sidebands.  However,  the  instability  should  be  treated  using  a  fully 
spatially  dependent  analysis.  Rather  than  using  an  initial  value  treatment  of 


>  t  -■ 


the  instability,  a  boundary  value  approach  would  be  more  appropriate,  as  the 
instability  is  convective  rather  than  absolute.  One  should  note  that  an  exact 
comparison  between  our  theory  and  the  experimental  results  of  Starke  and 
Malaberg  [1978]  is  only  valid  in  this  limit.  However,  even  with  this  simple 
approach,  one  can  see  the  similarities  between  the  broad  sidebands  they 
measure  and  the  ones  we  calculate.  One  could  also  note  a  similarity  between 
the  sidebands  of  the  present  analysis  and  the  "broad  bump"  in  the  plasma  line 
spectrum  observed  by  Thomson  radar  returns  [Showen  and  Kim,  1978]  during 
ionospheric  heating  experiments.  In  this  experiment,  the  wavelength  of  the 
Thomson  radar  is  fixed,  therefore,  the  broad  bump  is  made  of  only  those 
sidebands  that  match  with  the  radar  wavenumber.  The  recent  rocket  experiment 
of  Rose,  et  al.  [1985],  also  shows  evidence  of  broad  sideband  Langmuir  waves 
around  the  heater  frequency  Wo  measured  on  the  upleg  of  the  rocket  trajectory, 
below  the  reflection  layer. 
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3.  DISCUSSION 


The  second  order  analysis  presented  in  this  work  has  shown  that  the 
nonlinear  modification  of  the  background  distribution  function  by  a  resonant 
electrostatic  field  may  cause  tail  heating,  bump  production,  and  broad 
sideband  emissions.  The  power  levels  necessary  to  trigger  these  effects  in 
the  ionosphere's  &  and  F  regions  are  illustrated  by  the  following  example. 
Assuming  the  same  background  conditions  as  for  Figure  4,  i.e.,  A  -  0.01, 

Vf  «  10ve,  a  *  20  ve,  ve  ■  0.15  eV,  and  T  ~  10“ 5.  The  resulting  distribution 
is  plotted  in  Figure  5a  for  three  cases;  namely,  E0  -  0,  3ES,  and  6ES,  where 
Es  -  3.39/Ln  (in  SI  units)  for  the  preceeding  parameters.  The  threshold 
conditions  for  both  positive  slope  [Equation  (24)]  and  growth  [Equation  (28)] 
are  plotted  in  Figure  5b  and  appear  as  one  curve  since  the  two  plots  are 
almost  superimposed,  showing  that  for  F  region  parameters,  the  threshold  for 
growth  is  not  much  larger  than  that  for  positive  slope.  Outside  the  region 
contained  within  the  dotted  lines,  the  condition  of  (25)  is  not  satisfied  and 
the  denominators  of  both  (24)  and  (28)  become  negative  precluding  the 
formation  of  a  positive  slope  in  these  regions.  The  minimum  point  on  the 
curve  of  5b  is  determined  numerically  because  the  functions  involved  are 
transcendental.  In  this  example,  one  finds  E0,thr  “  4.83E3.  Therefore,  for 
the  plotted  value  of  E0  *  6ES,  one  is  well  above  threshold,  and  a  bump  does 
develop.  The  growth  rate  corresponding  to  this  value  of  E0  and  obtained  from 
(27)  is  plotted  in  figure  5c,  where  it  is  clear  that  growth  (i.e.,  y  >  0) 
occurs  only  in  the  region  where  the  distribution  has  a  positive  slope,  and 
that  it  maximizes  at  the  velocity  where  the  slope  is  maximum  in  the 


distribution.  Velocities  where  y  <  0,  correspond  to  damped  waves.  Assuming 


Mq  to  be  of  the  order  of  few  MHz,  which  is  usually  the  case  for  overdense 
heating  of  both  E  and  F  regions  in  the  ionosphere,  the  rise  time  for  the 


instability  is  only  a  few  hundred  microseconds.  Frequencies  at  the  edges  of 
the  sideband  grow  at  slower  rates.  In  Figure  5d,  curve  £  shows  the  width  Awg 
of  the  sidebands  spectrum  produced  for  a  homogeneous  medium  using  (29),  z  «  0. 
One  may  extend  the  validity  of  curve  £  to  an  inhomogeneous  medium  with  a  very 
large  density  scalelength  such  that  z/Ln  <<  1.  Curve  £  shows  the  other  limit 
case  wherein  z/Ln  ■  0.1,  a  finite  value  which  makes  the  condition  z  <  L(u) 
just  valid  for  the  ionospheric  plasma.  It  should  be  kept  in  mind  that  at  each 
value  of  z,  such  that  Az^  «  z  <  L(u),  the  distribution  function  with  the  bump 
exists  and  sidebands  are  being  produced.  This  means  that  curves  £  and  £  in 
Figure  5d  should  not  be  viewed  as  two  discrete  cases,  but  rather  as  two  bounds 
for  the  possible  spread  of  sidebands  ( 6013)  that  this  particular  example  would 
generate. 

To  justify  the  value  of  z  used  to  generate  curve  £  of  Figure  5d  in  terms 
of  slowing  down  length,  we  mention  that  for  energies  U(u)  of  the  electrons 
being  considered,  which  vary  usually  form  a  few  to  a  few  hundred  electron 
volts,  the  slowing  down  length  L(u),  evaluated  using  the  continuous  slowing 
down  approximation  where  L(u)  ■  U/(-  dU/dz),  has  a  minimum  value  that  varies 
from  a  few  kilometers  near  the  ionosphere's  upper  E  region  to  over  a  few 
hundred  kilometers  in  the  F  region.  This  result  is  displayed  in  Figure  6 
which  Incorporates  all  the  Important  slowing  down  processes  present  in  the 
ionosphere's  upper  E  and  F  regions.  Density  scale lengths ,  on  the  other  hand, 
vary  from  a  few  kilometers  in  the  E  region  to  one  hundred  kilometer  in  the  F 
region.  Therefore,  a  value  of  z  ~  0.1  Ln  <  L(u)  is  reasonable  for  an  upper 
range  of  validity  of  the  distribution  with  the  bump. 

Since  the  ionosphere  is  a  highly  variable  medium,  we  present  in  Table  1 
the  different  scalelength  dependent  parameters.  For  the  resonant  field's 


effective  phase  velocity  we  assume,  a  -  \>eLn  *  20ve,  there  exists  different 
combinations  for  the  product  Since  both  the  ionospheric  density 

scalelength  and  the  bulk  electron  collision  rate  depend  upon  altitude  and 
local  time,  we  compute  the  threshold  parameters  for  different  combinations  of 
vfi  and  Ln.  The  quantities  calculated  in  Table  1  use  the  same  parameters  that 
are  used  in  generating  the  previous  figures.  Typical  values  of  ve  and  Ln  in 
the  ionosphere's  upper  E  and  F  regions  have  a  product  that  varies  around  the 
value  chosen  ln  the  present  example. 

For  a  given  threshold  field  (E<j,thr)*  there  exists  a  corresponding 
energy  flux  density  for  the  electrostatic  wave  Pes  that  is  related  to  power 
delivered  by  the  high  frequency  electromagnetic  waves  at  the  reflection  layer 
pgp.  In  Table  1,  Pgp  is  estimated  assuming  a  mode  conversion  efficiency, 

1  *  pes/pHF  *  15%.  That  is,  Pyp  is  the  necessary  power  density  needed  at  the 
reflection  layer  (usually  anywhere  between  100  Km  and  350  Km  altitude)  in 
order  to  produce  an  electrostatic  field  just  above  the  threshold  value.  The 
values  displayed  in  Table  1  are  to  be  compared  with  vacuum  power  densities  of 
0.50  mW/m2  delivered  at  250  Km  altitude  by  the  Tromso  Heating  facility 
[Stubbe,  et  al.,  1982],  0.15  mW/m2  by  the  Alaska  HIPAS  facility  [Wong,  et  al., 
1983],  and  0.05  mW/m2  by  the  Arecibo  heater  [Mantas,  et  al.,  1981].  For  the 
example  we  are  presenting,  as  well  as  most  other  cases,  the  thresholds  for 
instability  are  within  these  values. 

It  is  also  of  interest  to  compare  the  thresholds  for  bump  instability 
that  we  have  calculated  with  those  necessary  to  produce  parametric  instabili¬ 
ties,  another  second  order  nonlinear  process  commonly  encountered  during 
ionospheric  modification  experiments.  The  two  mechanisms  pertinent  to  this 


situation  would  be  the  decay  instability  with  the  electromagnetic  pump 
decaying  near  its  cutoff  into  a  plasma  wave  and  an  ion  acoustic  wave,  and  the 
oscillating  two  stream  instability  involving  a  purely  growing  ion  mode.  The 
threshold  estimated  for  the  decay  instability  in  the  ionosphere  [Fejer,  1979] 
corresponds  to  an  HF  pump  amplitude  of  0.25  V/m  or  a  vacuum  HF  power 
PHF  *  180uW/m2  for  F  region  parameters.  For  the  oscillating  two  stream 
instability,  the  threshold  is  even  higher.  This  implies  that  the  sideband 
instability  thresholds  calculated  are  less  than,  and  at  most  on  the  order  of, 
the  smallest  parametric  instability  threshold,  namely  the  decay  process. 
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CONCLUSIONS 


The  second  order  analysis  presented  has  shown  the  following  major  effects 
in  response  to  the  excitation  of  a  spatially  localized  electric  field 
parallel  to  the  earth's  magnetic  field. 

1)  For  the  two  field  structures  "typical"  of  a  plasma  resonance,  the 
wavenumber  spectrum  is  broadband  and  unidirectional.  The  power 
spectra  of  both  the  Airy  and  the  Lorentzian  structures  are  identi¬ 
cal,  i.e.,  in  the  second  order  time  averaged  limit,  the  asymptotic 
modifications  in  the  distribution  function  are  identical.  However, 
one  should  note  that  the  Airy  pattern  has,  in  addition,  oscillating 
structures  in  it.  Its  equivalence  to  the  Lorentzian  is  therefore 
limited  to  situations  where  its  damped  envelope  is  the  main 
contribution  to  the  interaction,  which  is  the  case  in  the  present 
second  order  analysis.  Under  different  conditions,  the  oscillating 
part  of  the  Airy  structure  could  be  responsible  for  other  types  of 
interactions  not  considered  here  (e.g.,  particle  trapping). 

2)  Particle  acceleration  caused  by  direct  plasma  resonance  electric 
fields  (with  unidirectional  k-spectrum)  occurs  only  for  particles 
impinging  from  high  altitude  as  is  clearly  evident  from  (21). 

No  up  going  particles  are  resonantly  accelerated  by  such  a  field. 
This  result  implies  that  any  experimental  evidence  for  fast 
particle  acceleration  upstream  from  the  resonance  layer  indicates 
either  strong  pitch  angle  scattering  or  the  existence  of  electric 
fields  not  produced  by  direct  plasma  resonance  but  rather  by  some 
other  mechanism,  like  parametric  instabilities. 
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3)  The  Call  of  the  electron  distribution  is  heated  producing  an 
enhanced  fast  electron  flux  in  the  range  from  a  few  eV  to  few 
hundred  eV.  The  low  energy  part  of  this  tail  could  be  responsible 
for  the  observation  of  enhanced  airglow  at  63001  and  557 7  A 
[Sipler  and  Blondi,  1972],  and  may  account  for  the  high  energy 
particles  that  are  also  observed  [Carleson,  Wickwar,  and  Mantas, 
1982]. 

4)  Present  day  ionospheric  heating  experiments  appear  to  deliver 
powers  that  are  large  enough  to  overcome  the  threshold  for  bump 
formation.  Therefore,  a  bump  is  expected  to  develop  in  the  tail  of 
the  distribution  function.  Also,  it  is  found  that  for  the 
ionosphere,  this  threshold  is  quite  close  to  the  wave  instability 
threshold,  so  that  enhanced  noise  levels  should  also  occur. 

5)  Broad  sidebands  are  generated  in  the  F  region  by  the  bump  instabi¬ 
lity,  their  growth  time  is  of  the  order  of  10“ 1  millisecond.  They 
are  produced  at  different  heights  downstream  from  the  interaction 
region,  up  to  the  slowing  down  length  or  nonlinear  saturation 
length,  whichever  is  shortest.  For  a  distant  observer,  the  net 
broad  sideband  spectrum  consists  of  the  integrated  effect  from  all 
these  heights.  Experiments  measuring  the  plasma  line  spectrum 
using  Thomson  radar  backscattering  from  the  ionosphere  during  HF 
heating  {Kantor,  1974;  Showen  and  Kim,  1978],  show  evidence,  under 
some  conditions,  of  a  "broad  bump”  in  the  spectrum.  In  addition, 
the  broad  spectrum  of  plasma  waves  around  the  HF  heater  frequency 
measured  in  the  HERO  rocket  flights  [Rose,  et  al.,  1985],  could  be 
another  example  of  the  phenomena  considered  by  the  present  study. 


In  che  laboratory,  Che  results  of  Starke  and  Malmberg  [1978]  have 
simularlties  to  ours,  both  In  the  broad  sideband  spectrum  they  find 
and  In  the  quaslllnear  theory  they  Invoke.  However,  since  the 
spatial  dependence  of  che  growing  sidebands  Is  not  self-consistently 
Included  In  our  analysis,  a  direct  comparison  Is  not  possible  at 
this  point. 


APPENDIX:  FOURIER  TRANSFORMS  OF  THE  AIRY  PATTERN  AND  THE  LORENTZIAN 


We  derive  here  equaClons  (17)  and  (18).  The  Fourier  transform  is  defined 


E(k)  ■  /  dz  Ej(z)  exp(-ikz) 


where  Ej(z)  is  defined  by  (15)  for  the  Airy  pattern  and  (16)  for  the 
Lorentzian. 


One  may  write  (15)  as  [Abramowltz  and  Stegun,  1965] 

E|(z)  •  E0(kDLn)2/3/  dt  exp  {-1  [t-  +  (kDLn)2/3  (-  f - IT)  t  +  t]}.(A2) 

o  n  *- 

Substituting  in  (Al)  and  integrating  over  all  space 

EA(k)  -  2irE0(kDLn)2/3  dt  exp  {-  i[^ - (koLr)2/3rt  +  y]}«5[k  -  at] 

*  2irEoLn  exp[-k/kQ  -  i£-y  +  f)],  k  >  o  (A3) 

ja  *• 


,  k  <  o 


where  k0  -  (n^)-!,  a  *  (koLn)2/3/Ln»  and  kg  Is  the  Deybe  wavenumber. 
Equations  (A3)  and  (A4)  are  identical  to  (17). 


For  the  Lorentzian  case,  (Al)  is  rewritten  as 


rL(k)  -  /  < dz  - -rnr  exp(*ikz) 


2*  E0Ln  exp(-k/kQ  -  iir/2 ) ,  k  >  0 


toy 

ran 

I 


| 


m 
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where  Che  integration  is  carried  out  in  the  complex  z  plane  closing  the 
contour  in  Che  lower  half-plane  for  k  >  0  (encircling  the  pole  at  -in^),  and 
in  the  upper  half-plane  for  k  <  0  (where  no  poles  are  enclosed).  Again, 
equations  (A5)  and  (A6)  are  identical  to  (18). 
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FIGURE  CAPTIONS 

Figure  1.  Schematic  diagram  Illustrating  the  geometry  of  the  problem. 

Electrons  from  the  high  density  side  (on  the  left)  pass  through 
the  localized  interaction  region  at  plasma  resonance  (z  ■  0), 
having  a  characteristic  width  Az^.  They  emerge  accelerated  in  the 
region  zQ  £  z  «  L(u)  forming  an  enhanced  tail  and,  above  a 
threshold  level,  a  bump  in  the  distribution.  The  bump  eventually 
relaxes  to  the  equilibrium  background  distribution  (on  the  far 
right)  after  a  distance  longer  than  a  slowing  down  length  L(u). 

Figure  2.  Spatial  dependence  of  resonant  parallel  electric  fields  for  the 

Airy  (continuous  lines)  and  the  Lorentzian  (dashed  lines)  patterns 
for  different  damping  values:  (a)  Magnitudes  for  F'  -  0.2  and  0.8; 
(b)  Real  parts  for  T*  -  0.4.  Note  the  wavelike  structure  of  the 
Lorentzian  in  the  neighborhood  of  the  origin. 

Figure  3.  Dual  plot  of  the  magnitude  of  the  Fourier  transformed  field  (for 
either  the  Airy  or  the  Lorentzian  pattern)  versus  wavenumber  k 
(upper  scale  and  dashed  curve),  and  versus  velocity  u  *  o>0/k  (lower 
scale  and  continuous  curve).  Normalization  is  such  that 
kQ  ■  w0/a  «  (rLjj)-1,  a  is  the  effective  phase  velocity  of  the 
resonant  field. 


Figure  4.  Modification  in  electron  velocity  distribution  obtained  from 
equation  (22)  for  E0  *  3Ea,  a  -  2vf  -  20ve,  and  A  -  It. 
Unperturbed  background  distribution  (E0  ■  0)  is  plotted  for 
comparison. 


Figure  5.  Production  of  a  bump  in  the  distribution  function  and  associated 
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Figure  6 


sideband  instability:  (a)  the  distribution  function  with  the  tail 
magnified  (scale  to  the  right)  for  a  ■  2vj  »  20ve,  ve  ~  0.15  eV, 

A  ■  IX,  a>0/2x  ~  5  MHz,  r  -  10"5,  and  different  values  of  E0;  (b) 
threshold  curves  for  bump  formation  and  instability  generation 
(both  curves  superimposed);  (c)  temporal  growth  rates  for  the 
unstable  sidebands  (y  >  0)  generated  on  the  positive  slope  region 
of  the  bump  for  E0  »  6ES;  (d)  frequency  range  Ao»g  corresponding  to 
sidebands  generated  locally  in  a  slowly  varying  medium  for  curve 
a.:  z  m  0,  and  curve  J>:  z  »  0.1Ln.  A  distant  observer  downstream 
from  the  resonance  region  may  sample  sidebands  within  a  frequency 
range  5uig. 

Energy  dependence  of  slowing  down  length  L(u)  for  fast  electrons  in 
the  upper  E  region  (at  150  Km  altitude)  and  the  F  region  (at  300 
Km  altitude).  Estimates  are  based  on  a  model  incorporating  the 
most  dominant  slowing  down  mechanisms  of  the  ionosphere:  impact 
excitations,  ionizations,  recombinations,  electron-electron 
collisions,  •  .  .  [see,  for  example,  Stammes  and  Rees  (1983)  and 
the  references  therein] .  The  background  ionosphere  used  is  the  one 
given  by  Shoucri,  Morales,  and  Maggs  [1984]. 


Energy  (eV) 
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TABLE  1.  Electric  Field  and  Power  Densities  Thresholds  for  Different 
Ionospheric  Parameters* 


Ln(Km) 

5 

10 

20 

50 

Vs"1) 

920 

460 

230 

92 

Eo,thr(“v*“'1) 

3.27 

1.64 

0.82 

0.33 

Pe8(uW.nr2) 

23.6 

11.8 

5.89 

2.36 

pHp( uW.m-2) 

157 

78.7 

39.3 

15.7 

*For  a  »  20ve  ■  vgLn,  various  possible  combinations  of  ionospheric 
density  scalelengths  L„,  and  bulk  electron  collisions  vg  yield: 
threshold  field  of  (23)  for  sideband  instability  E0>thr;  threshold 
energy  flux  density  in  the  electrostatic  plasma  wave  Pes;  and 
corresponding  threshold  vacuum  HF  power  at  the  reflection  layer  P^p 
assuming  a  mode  conversion  efficiency  of  n  *  pes/pHF  "  °»15. 
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